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Abstract 

Ferromagnetic resonance in thin films is analyzed under the influence of spatiotemporal feedback 
effects. The equation of motion for the magnetization dynamics is nonlocal in both space and time 
and includes isotropic, anisotropic and dipolar energy contributions as well as the conserved Gilbert- 
and the non-conserved Bloch-damping. We derive an analytical expression for the peak-to-peak 
linewidth. It consists of four separate parts originated by Gilbert damping, Bloch-damping, a mixed 
Gilbert-Bloch component and a contribution arising from retardation. In an intermediate frequency 
regime the results are comparable with the commonly used Landau-Lifshitz-Gilbert theory combined 
with two-magnon processes. Retardation effects together with Gilbert damping lead to a linewidth 
the frequency dependence of which becomes strongly nonlinear. The relevance and the applicability 
of our approach to ferromagnetic resonance experiments is discussed. 
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I. INTRODUCTION 



Ferromagnetic resonance enables the investigation of spin wave damping in thin or ul- 
trathin ferromagnetic films. The relevant information is contained in the linewidth of the 
resonance signal [T]-[3]. Whereas the intrinsic damping included in the Gilbert or Landau- 
Lifshitz-Gilbert equation [H |5], respectively, predicts a linear frequency dependence of the 
linewidth [6j, the extrinsic contributions associated with two-magnon scattering processes 
show a nonlinear behavior. Theoretically two-magnon scattering was analyzed for the case 
that the static external field lies in the film plane [71 [8]. The theory was quantitatively 
validated by experimental investigations with regard to the film thickness [9j. Later the 
approach was extended to the case of arbitrary angles between the external field and the 
film surface [lOj. The angular dependence of the linewidth is often modeled by a sum of 
contributions including angular spreads and internal field inhomogeneities [TT]. Among oth- 
ers, two-magnon mechanisms were used to explain the experimental observations |T2l - [T7] 
whereas the influence of the size of the inhomogeneity was studied in |18) . As discussed in 
[21 [13] the two-magnon contribution to the linewidth disappears for tipping angles between 
magnetization and film plane exceeding a critical one = 7r/4. Recently, deviations from 
this condition were observed comparing experimental data and numerical simulations |I7] . 
Spin pumping can also contribute to the linewidth as studied theoretically in [19j. How- 
ever, a superposition of both the Gilbert damping and the two-magnon contribution turned 
out to be in agreement very well with experimental data illustrating the dependence of the 
linewidth on the frequency [T6l I20H23] . Based on these findings it was put into question 
whether the Landau-Lifshitz-Gilbert equation is an appropriate description for ferromag- 
netic thin films. The pure Gilbert damping is not able to explain the nonlinear frequency 
dependence of the linewidth when two-magnon scattering processes are operative [31 
Assuming that damping mechanisms can also lead to a non-conserved spin length a way 
out might be the inclusion of the Bloch equations [25] [26] or the the Landau-Lifshitz-Bloch 
equation [27[ [28] into the concept of ferromagnetic resonance. 

Another aspect is the recent observation [29] that a periodic scattering potential can alter 
the frequency dependence of the linewidth. The experimental results are not in agreement 
with those based upon a combination of Gilbert damping and two-magnon scattering. It 
was found that the linewidth as function of the frequency exhibits a non monotonous be- 
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havior. The authors |29) suggest to reconsider the approach with regard to spin relaxations. 
Moreover, it would be an advantage to derive an expression for the linewidth as a measure 
for spin damping solely from the equation of motion for the magnetization. 

Taking all those arguments into account it is the aim of this paper to propose a gener- 
alized equation of motion for the magnetization dynamics including both Gilbert damping 
and Bloch terms. The dynamical model allows immediately to get the magnetic susceptibil- 
ity as well as the ferromagnetic resonance linewidth which are appropriate for the analysis 
of experimental observations. A further generalization is the implementation of nonlocal 
effects in both space and time. This is achieved by introducing a retardation kernel which 
takes into account temporal retardation within a characteristic time r and a spatial one 
with a characteristic scale ^. The last one simulates an additional mutual interaction of 
the magnetic moments in different areas of the film within the retardation length ^. Re- 
cently such nonlocal effects were discussed in a complete different context j30]. Notice that 
retardation effects were already investigated for simpler models by means of the Landau- 
Lifshitz-Gilbert equation. Here the existence of spin wave solutions were in the focus of the 
consideration ^31j. The expressions obtained for the frequency/damping parameters were 
converted into linewidths according to the Gilbert contribution which is a linear function 
of the frequency fST\ |32]- In the present approach we follow another line. The propagating 
part of the varying magnetization is supplemented by the two damping terms due to Gilbert 
and Bloch, compare Eq. Based on this equation we derive analytical expressions for the 
magnetic susceptibility, the resonance condition and the ferromagnetic resonance linewidth. 
Due to the superposition of damping and retardation effects the linewidth exhibits a non- 
linear behavior as function of the frequency. The model is also extended by considering 
the general case of arbitrary angles between the static external field and the film surface. 
Moreover the model includes several energy contributions as Zeeman and exchange energy 
as well as anisotropy and dipolar interaction. The consequences for ferromagnetic resonance 
experiments are discussed. 

II. DERIVATION OF THE EQUATION OF MOTION 

In order to define the geometry considered in the following we adopt the idea presented 
in [To], i.e. we employ two coordinate systems, the xyz-sjstem referring to the film surface 
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FIG. 1. (Color online) The geometry referring to the film and the magnetization. Further descrip- 
tion in the text. 

and the XYZ-system which is canted by an angle Bm with respect to the film plane. The 
situation for a film of thickness d is sketched in Fig. [1} The angle 0m describing the direction 
of the saturation magnetization, aligned with the Z-axis, originates from the static external 
field Ho which impinges upon the film surface under an angle Gh- Therefore, it is more 
convenient to use the XYZ-system for the magnetization dynamics. As excitation source 
we consider the radio-frequency (rf) magnetic field h^f pointing into the x = X-direction. It 
should fulfill the condition hj-f ^ ifo- To get the evolution equation of the magnetization 
M(r, t), r = {x, y, z) we have to define the energy of the system. This issue is well described 
in Ref. [TU], so we just quote the most important results given there and refer to the cited 
literature for details. Since we consider the thin film limit one can perform the average along 
the direction perpendicular to the film, i.e. 

1 rd/2 

M(r||,t) = - / d|/M(r,t), (1) 

" J-d/2 

where ry = (x, 0, z) lies in the film plane. In other words the spatial variation of the 
magnetization across the film thickness d is neglected. The components of the magnetization 
point into the directions of the XYZ-system and can be written as [33] 

M(r||,t) = Mx(r||)ex + MY(r||)eY+ ( Ms- "^^^ "M ez . (2) 
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Typically the transverse components Mx,y are assumed to be much smaller than the satu- 
ration magnetization Ms. Remark that terms quadratic in Mx,y in the energy will lead to 
linear terms in the equation of motion. The total energy of the system can now be expressed 
in terms of the averaged magnetization from Eq. ([T]) and reads 

"H = "Hz + "Hex + "Ha + "Hd . (3) 
The different contributions are the Zeeman energy 

Hz = - y"dVi7osin(eH-eM)MY(r||) 



the exchange energy 

"^^^^ 2Ms 



J dV [VMx(rij)]' + [VMY(r|,)]' , (5) 



the surface anisotropy energy 

Ha =;^^ sin^ (Gm) + ^ sin(2eM) J d'rMy{r^^) 

+ ^ cos(2eM) j dVMY(r|i)2-sin2(eM) j dVMx(r|| 



and the dipolar energy 



«.^2.M|V.„=(e„) + ./d3.{2M.s,.(2eM)Mv(r,) 

+ i^-^ sm^iBu) - {dk\\ - 2) cos2(eM) - 2 sin2(eM) j MY(r|| 
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+ ( ^ - 2 sin2(eM)) Mx(r||)2 - sin(eM) Mx(r|,)MY(r|| 



(6) 



(7) 



In these expressions V = Ixhd is the volume of the film, D designates the exchange stiffness 
and ifs oc represents the uniaxial out-of-plane anisotropy field. If < the easy axis 
is perpendicular to the film surface. The in-plane anisotropy contribution to the energy is 
neglected but it should be appropriate for polycrystalline samples jT6]. Moreover = |k|j| 
is introduced where ky = k^e^ + kz is the wave vector of the spin waves parallel to the 
film surface. Eqs. (|3])-([7]) are valid in the thin film limit k\\d <^ 1. In order to derive "Hd in 
Eq. ([T]) one defines a scalar magnetic potential and has to solve the corresponding boundary 



value problem inside and outside of the film |3l]. As result [10] one gets the expressions in 
Eq. 0. 

In general if the static magnetic field is applied under an arbitrary angle Bh the mag- 
netization does not align in parallel, i.e. Bm 7^ Bh. The angle Bm can be derived from 
the equilibrium energy "Heq = 'H{Mx = 0,My = 0). Defining the equilibrium free energy 
density as fcqiQm) = T^eq/V according to Eqs. (|3])-([7]) one finds the well-known condition 

sin(BH - Bm) = ^^^^ + ^^^^2 Bm) (8) 

by minimizing /cq with respect to Bm. We further note that all terms linear in My in 
Eqs. (|3|-([7]) cancel mutually by applying Eq. Q as already pointed out in Ref. [10]. 

The energy contributions in Eqs. ([s]) and the geometric aspects determine the dynamical 
equation for the magnetization. The following generalized form is proposed 



^M(r,|,t)= // drj|dt'r(r||-rj|;t-0<!7 [Heff(rj|, t') x M(rj|, t' 



dt 



+ a 



M(rj|,t')x|;M(rj|,tO 
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(9) 



where 7 = gfis/^ is the absolute value of the gyromagnetic ratio, T2 is the transverse 
relaxation time of the components M_l = Mx ex + My ey and a denotes the dimensionless 
Gilbert damping parameter. The latter is often transformed into G = a'yMs representing the 
corresponding damping constant in unit s~^. The effective magnetic field Hes is related to 
the energy in Eqs. ([3])-([7]) by means of variational principles [35|, i.e. Heg = — 5'H/(5M + hi.f. 
Here the external rf-field h^f{t) is added which drives the system out of equilibrium. 

Regarding the equation of motion presented in Eq. (|9| we note that a similar type was 
applied in [12] for the evaluation of ferromagnetic resonance experiments. In this paper 
the authors made use of a superposition of the Landau-Lifshitz equation and Bloch-like 
relaxation. Here we have chosen the part which conserves the spin length in the Gilbert form 
and added the non-conserving Bloch term in the same manner. That the combination of 
these two distinct damping mechanisms is suitable for the investigation of ultrathin magnetic 
films was also suggested in [2^. Since the projection of the magnetization onto the Z-axis is 
not affected by T2 this relaxation time characterizes the transfer of energy into the transverse 
components of the magnetization. This damping type is supposed to account for spin-spin 
relaxation processes such as magnon-magnon scattering [33| 136] . In our ansatz we introduce 



another possible source of damping by means of the feedback kernel r(r|| — rj|; t — t'). The 
introduction of this quantity reflects the assumption that the magnetization M(r||,t2) is 
not independent of its previous value lSA{r\\,ti) provided t2 — ti < r. Here r is a time 
scale where the temporal memory is relevant. In the same manner the spatial feedback 
controls the magnetization dynamics significantly on a characteristic length scale ^, called 
retardation length. Physically, it seems to be reasonable that the retardation length differs 
noticeably from zero only in 2;-direction which is shown in Fig. [TJ As illustrated in the figure 
Nl{x, Zi,t) is affected by l\/l{x, Z2,t) while M(x, 2:3,^) is thought to have negligible influence 
on lSA{x, Zi,t) since {z^ — Zi\ > ^. Therefore we choose the following combination of a local 
and a nonlocal part as feedback kernel 



rfni 



r'ft 



t') 



ro5(r|| 
To 



+ 



4^r 



-rj|)5(t-t') 
6{x — x') exp 



~\z 



exp 



it-t') 



t > t' 



(10) 
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The intensity of the spatiotemporal feedback is controlled by the dimensionless retardation 



strength Tq. The explicit form in Eq. (10) is chosen in such a manner that the Fourier- 
transform r(k||,w) — )■ Fq for ^ — )■ and r — )■ 0, and in case Fq = 1 the ordinary equation 
of motion for the magnetization is recovered. Further, J dr||dtF(r||,t) = Fq < 00, i.e. the 
integral remains finite. 



III. SUSCEPTIBILITY AND FMR-LINEWIDTH 

If the rf-driving field, likewise averaged over the film thickness, is applied in X-direction, 
i.e. hi.f(r||,t) = hx{'r\\,t) ex, the Fourier transform of Eq. ^ is written as 



_7F(k||,u;) 7T2 



i^2i(kj| 



7F(k||,w) 7T2 



+ ^ + ^i2(k|| 



Mx(k||,u;) 



MY(k|,,u;) 



^i(k|| 



^2(k||) + 
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MY(k||,w) 
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Mx(k||,a;) -Ms/ix(k||,a;). 

(11) 
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The effective magnetic fields are expressed by 

i/i(k||) =Ho cos(eH - Om) + (47rMs + Hs) cos(2 Gm) 

+ 27idk\\Ms sin2(eM) - cos2(eM)^ + D kl 

i72(k||) =Ho cos(eH - Om) - (47rMs + Hs) sin^Qu) 

+ 2TidMs -r^ + Dk^, 

k\\ " 



(12) 



and 



i/i2(k|j) = 2'ndMs ^ sin(eM) = -i^2i(k||) . (13) 



The Fourier transform of the kernel yields 

Ti To + Ti i ^ 

— ~ Tiwr, 

(14) 

where the factor 1/2 arises from the condition t > t' when performing the Fourier trans- 



^ , Fo (1 + iut) + Fi Fo + Fi i 

2(l + ic.r) " ^ 2^^""^ 



formation from time into frequency domain. In Eq. (14) we discarded terms uj t <^ 1. 
This condition is fulfilled in experimental realizations. So, it will be turned out later the 
retardation time r ~ 10 fs. Because the ferromagnetic resonance frequencies are of the order 
10 . . . 100 GHz one finds w^r^ ~ 10~^...10^^. The retardation parameter P = ^k^, introduced 



in Eq. (14), will be of importance in analyzing the linewidth of the resonance signal. With 



regard to the denominator in Fi, compare Eq. (14), the parameter (3 may evolve ponderable 
influence on the spin wave damping if this quantity cannot be neglected compared to 1. 
As known from two-magnon scattering the spin wave modes can be degenerated with the 
uniform resonance mode possessing wave vectors k\\ ~ lO^cm"^. The retardation length ^ 
may be estimated by the size of inhomogeneities or the distance of defects on the film sur- 
face, respectively. Both length scales can be of the order ~ 10...1000nm, see Refs. [T5| 
Consequently the retardation parameter /3 could reach or maybe even exceed the order of 1. 

Let us stress that in case /3 = 0, r = 0, Fq = 1 and neglecting the Gilbert damping, 
i.e. a = 0, the spin wave dispersion relation is simply Hi(k\^)H2{'k\\) — i/i2(k||)- This 
expression coincides with those ones given in Refs. |7| and [TOj . 

Proceeding the analysis of Eq. ([TT| by defining the magnetic susceptibility x as 



Ma{kii,u) = ^Xaf3{kii,u)h(s{ki\,uj) , = {X, Y} , (^^5) 

/3 
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where h/^ plays the role of a small perturbation and the susceptibility Xa/s exhibits the 



response of the system. Eq. ( 15 ) reflects that there appears no dependence on the direction 
of k||. 

Since the rf-driving field is applied along the ex-direction it is sufficient to focus the 



following discussion to the element Xxx ^^e susceptibility tensor. From Eq. (11) we 
conclude 



Xxx(^||, 



Ms 



Hi{kii,uj) + 



7 



7 



Ho 



+ 



(16) 



Because at ferromagnetic resonance a uniform mode is excited let us set k\\ = in Eqs. (12)- 



(13). Considering the resonance condition we can assume /3 = ^/c^ = 0. For reasons men- 
tioned above we have to take /3 = C,kz ^ when the linewidth as a measure for spin damping 
is investigated. Physically we suppose that spin waves with non zero waves vectors are not 
excited at the moment of the ferromagnetic resonance. However such excitations will evolve 



during the relaxation process. In finding the resonance condition from Eq. (16) it seems to 
be a reasonable approximation to disregard terms including the retardation time r. Such 
terms give rise to higher order corrections. In the same manner all the contributions orig- 
inated from the damping, characterized by a and T2, are negligible. Let us justify those 
approximation by quantitative estimations. The fields Hi, H2 and u/'~f are supposed to 
range in a comparable order of magnitude. On the other hand one finds a ~ 10^'^. ..10^^, 
UT2 ~ 10~^ and UT ~ 10~^. Under these approximations the resonance condition reads 



(^^^ =rlHi{k\\=0)H2{k^\=0). 



(17) 



This result is well known for the case without retardation with Fq = 1. Although the retarda- 
tion time r and the retardation length ^ are not incorporated in the resonance condition, the 



strength of the feedback may be important as visible in Eq. (17). Now the consequences for 
the experimental realization will be discussed. To address this issue the resonance condition 



Eq. (17) is rewritten in terms of the resonance field H,. = Hq{uj = Uj.) leading to 

1 



Hr 



2 cos(0H - 




(47rMs + Hs)^ cos4(0m) + 



1 2w/ 
ro 7 . 



(18) 



-(47rMs + i7s)(l-3 sin2(eM)) 



9 



90 
80 
70 

„ 60 



§ 40 
® 30 
20 
10 


10 20 30 40 50 60 70 80 90 
Oh [deg] 

FIG. 2. (Color online) Dependence of the magnetization angle 0m on the angle Gh under which the 
static external field is applied for (x>r/ (27r) = 10 GHz. The parameters are taken from |16| : 47rMg = 
16980 G, Hs = -3400 G, 7 = 0.019 GHz/G. 




The result is arranged in the in the same manner as done in [16]. The difference is the 
occurrence of the parameter Fq in the denominator. In [16j the gyromagnetic ratio 7 and 
the sum (47rMg + Hg) were obtained from Bh -dependent measurements and a fit of the 



data according to Eq. (18) with Fq = 1 under the inclusion of Eq. (|8]). If the saturation 
magnetization can be obtained from other experiments [I6] the uniaxial anisotropy field Hg 
results. Thus, assuming Fq 7^ 1 the angular dependence 0m (©h) and the fitting parameters 
as well would change. In Fig. [2] we illustrate the angle 0m (©h) for different values of Fq and 
a fixed resonance frequency. If Fq < 1 the curve is shifted to larger Bm and for Fq > 1 to 
smaller magnetization angles. To produce Fig. [2] we utilized quantitative results presented 
in |16]. They found for Co films grown on GaAs the parameters AtcMs = 16980 G, Hs = 



—3400 G and 7 = 0.019 GHz/G. As next example we consider the influence of Hs and denote 
i7g°^ = —3400 G the anisotropy fleld for Fq = 1 and H^^^ the anisotropy field for Fq 7^ 1. The 
absolute value of their ratio \H^^^ /H^^\, derived from Hr{H^\ Fq = 1) = Hr{H^^\ Fq 7^ 1), 
is depicted in Fig. |3]for various frequencies. In this graph we assumed that all other quantities 
remain fixed. The effect of a varying retardation strength on the anisotropy field can clearly 
be seen. The change in the sign of the slope indicates that the anisotropy field H^ may even 
change its sign. From here we conclude that the directions of the easy axis and hard axis 
are interchanged. For the frequencies 4 GHz and 10 GHz this result is not observed in the 
range chosen for Fq. Moreover, the effects become more pronounced for higher frequencies. 
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FIG. 3. (Color onhne) Effect of varying retardation strength on the uniaxial anisotropy field for 
various frequencies and Qm = vr/3. AttMs = 16980 G, Hs = -3400 G, 7 = 0.019 GHz/G, see [l6]. 



In Fig. |3]we consider only a possible alteration of the anisotropy field. Other parameters like 
the experimentally obtained gyromagnetic ration were unaffected. In general this parameter 
may also experiences a quantitative change simultaneously with iYg. 



Let us proceed by analyzing the susceptibility obtained in Eq. (16). Because the following 



discussion is referred to the energy absorption in the film, we investigate the imaginary part 
of the susceptibility Xxx ■ Since experimentally often a Lorentzian curve describes sufficiently 
the resonance signal we intend to arrange Xxx i'^ ^^e form Aq/{1 + u^), where Aq is the 
absolute value of the amplitude and m is a small parameter around zero. The mapping to a 
Lorentzian is possible under some assumptions. Because the discussion is concentrated on 
the vicinity of the resonance we introduce 6H = Hq — Hj., where Hj. is the static external 



field when resonance occurs. Consequently, the fields in Eq. (12) have to be replaced by 
Hi^2 h[''1 + 6H cos(Oh — ©m)- Additionally, we take into account only terms of the order 
y/eX in the final result for the linewidth where {e. A} oc {uj/'y[a + lot] + 1/ (7T2)}. After a 
lengthy but straightforward calculation we get for SH/H^l <^ 1 and using the resonance 



condition in Eq. (17) 



1 + 



An 



At 



2 ' 



An 



H. 



(r) 



[I + k) COs(0H-eM)AT' h[ 



(r) 



(19) 



Here we have introduced the total half-width at half-maximum (HWHM) At which can be 
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brought in the form 



An 



A2 +A2 



A2 



(20) 



cos(eH - ©m) 

The HWHM is a superposition of the Gilbert contribution Aq, the Bloch contribution Ab, 
a joint contribution Aqb arising from the combination of the Gilbert and Bloch damping 
parts in the equation of motion and the contribution Ar which has its origin purely in the 
feedback mechanisms introduced into the system. The explicit expressions are 




a — 



16v^ Tori 

(i + /€) (ro + ri)3^ 



At 



4rn 



(ro + ri)(i + ft: 




4ri 



(ro + ri)2 7 7r2 



A 



Sir 



K 



GB 



A 



R 



(ro + ri)(i + /c) 72T2 
:i + /c)7(ro + ri)3- 



(21a) 
(21b) 

(21c) 
(21d) 



The parameter Fi is defined in Eq. (^14|. If the expressions under the roots in Eqs. (21a) 



and (21b) are negative we assume that the corresponding process is deactivated and does 



not contribute to the linewidth AHt- Typically, experiments are evaluated in terms of the 
peak-to-peak linewidth of the derivative dxxx/d-f^o, denoted as AHrj. One gets 

2 



AHr, 



(22) 



where the index 77 stands for G (Gilbert contribution), B (Bloch contribution), GB (joint 
Gilbert-Bloch contribution), R (pure retardation contribution) or T designating the total 



linewidth according to Eq. (20) and Eqs. (21a)-(21d). Obviously these equations reveal a 



strong nonlinear frequency dependence, which will be discussed in the subsequent section. 



IV. DISCUSSION 



As indicated in Eqs. (20) - (22) the quantity Aif^ consists of well separated distinct 
contributions. The behavior of Aif^ is shown in Figs. [4]- [6] as function of the three retardation 
parameters, the strength Fq, the spatial range /3 and the time scale r. In all figures the 
frequency / = uj/{2n) is used. In Fig. lilthe dependence on the retardation strength Fq is 
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FIG. 4. (Color online) Influence of the retardation strength Fq on the peak-to-peak linewidth Ai^T 
for various frequencies (top graph) and on the single contributions /S.H^ for / = 70 GHz (bottom 
graph). Ab = is this frequency region. The parameters are: Gh = 6m = 0, /3 = 0.5, a = 0.01, 
= 5 X 10-s s, T = 1.7 x 10^^"^ s. The other parameters are 47rMs = 16980 G, = -3400 G, 7 = 
0.019 GHz/G, compare P] . 

shown. As already observed in Figs. [2] and [3] a small change of Fq may lead to remarkable 
effects. Hence we vary this parameter in a moderate range 0.5 < Fq < 2. The peak-to-peak 
linewidth Ai^x as function of Fq remains nearly constant for / = 4 GHz and / = 10 GHz, 
whereas for / = 35 GHz a monotonous growth-up is observed. Increasing the frequency 
further to / = 50 GHz and 70 GHz the curves offers a pronounced kink. The subsequent 
enhancement is mainly due to the Gilbert damping. In the region of negative slope we 
set /\Hq{Y'q) = 0, while in that one with a positive slope AHq^Tq) > grows and tends 
to 2au/{^/3'y) for Fq — 00. The other significant contribution A/Jr, arising from the 
retardation decay, offers likewise a monotonous increase for growing values of the retardation 
parameter Fq. This behavior is depicted in Fig. |4]for / = 70 GHz. Now let us analyze the 
dependence on the dimensionless retardation length /3 = C,kz- Because /3 is only nonzero if 
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FIG. 5. (Color online) Influence of the dimensionless retardation length /3 = (^k^ on the total 
peak-to-peak linewidth AH'-£ for various frequencies (top graph) and on the single contributions 
AH^ for / = 70 GHz (bottom graph); Ab = in this range. The parameters are: Bh = Qm = 0, 
To = 1.1, a = 0.01, Ts = 5 X 10-*^s,t = 1.7 x 10"^^ s. The other parameters: AttMs = 16980 G, 
Hs = -3400 G and 7 = 0.019 GHz/G are taken from |I6]. 



^ this parameter ^ accounts the influence of excitations with nonzero wave vector. We 
argue that both nonzero wave vector excitations, those arising from two-magnon scattering 
and those originated from feedback mechanisms, may coincide. Based on the estimation 
in the previous section we consider the relevant interval 10^^ < < 10. The results are 
shown in FigjSj Within the range of /3 one recognizes that the total peak-to-peak linewidths 
AHt for / = 4 GHz and / = 10 GHz offer no alteration when /3 is changed. The plotted 
linewidths are characterized by a minimum followed by an increase which occurs when /3 
exceeds approximately 1. This behavior is the more accentuated the larger the frequencies 
are. The shape of the curve can be explained by considering the single contributions as 
is visible in the lower part in Fig. [s] While both quantities AHq{(3) and AH^((3) remain 
constant for small /3, AHq{/3) tends to a minimum and increases after that. The quantity 
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FIG. 6. (Color online) Influence of the retardation time r on the total peak-to-peak linewidth 
AHt for various frequencies (top graph) and on the single contributions Aff^ for / = 70 GHz 
(bottom graph). Ab = in this region. The parameters are ©h = ©m = 0, /? = 0.5, a = 0.01, 
T2 = 5 X 10~'^s,ro = 1.1; the other parameters are taken from [16]: AnMs = 16980 G, Hs = 
-3400 G, 7 = 0.019 GHz/G. 



A/7r(/3) develops a maximum around /3 ^ 1. Thus, both contributions show nearly opposite 
behavior. The impact of the characteristic feedback time r on the linewidth is illustrated 
in Fig. |6| In this figure a linear time scale is appropriate since there are no significant 
effects in the range Ifs > r > 0. The total linewidth AHt{t) is again nearly constant 
for / = 4 GHz and / = 10 GHz. In contrast Ai^xl^) reveals for higher frequencies two 
regions with differing behavior. The total linewidth decreases until AHq{t) becomes zero. 
After that one observes a positive linear slope which is due to the retardation part AH^ {t). 



This linear dependency is recognizable in Eq. (21d), too. Below we will present arguments 



why the feedback time r is supposed to be in the interval < r < 100 fs. Before let us 
study the frequency dependence of the linewidth in more detail. The general shape of the 
total linewidth AH-j^co) is depicted in Fig. [7| Here both the single contribution to the 



15 



350 
^ 300 

250 
<1 200 
^ 150 



20 40 60 80 100 120 
frequency / [GHz] 

FIG. 7. (Color online) Frequency dependence of all contributions to the peak-to-peak linewidth for 
Sh = Om = 0, /3 = 0.5, a = 0.01, T2 = 5 x 10"^ s, r = 1.7 x 10"^^ s and Fq = 1.2. Parameters taken 
from Ref. [E]: 47rMs = 16980 G, Hs = -3400 G and 7 = 0.019 GHz/G. The Bloch contribution 
AHb is shown in the inset. 




linewidth and the total linewidth are shown. Notice that the total linewidth is not simply 



the sum of the individual contributions but has to be calculated according to Eq. (20). One 
realizes that the Bloch contribution AH^ is only nonzero for frequencies f < Q GHz in the 
examples shown. Accordingly AHb = in Figs. |4]|6] (lower parts) since these plots refer to 
/ = 70 GHz. The behavior of the Gilbert contribution deviates strongly from the typically 
applied linear frequency dependence. Moreover, the Gilbert contribution will develop a 
maximum value and eventually it disappears at a certain frequency where the discriminant 



in Eq. (21a) becomes negative. Nevertheless, the total linewidth is a nearly monotonous 
increasing function of the frequency albeit, as mentioned before, for some combinations of 
the model parameters there might exist a very small frequency region where AHg reaches 
zero and the slope of AHt becomes slightly negative. The loss due to the declining Gilbert 
part is nearly compensated or overcompensated by the additional line broadening originated 
by the retardation part and the combined Gilbert-Bloch term. The latter one is AHq^ oc ^/f 



and AHji oc / , see Eqs. (21c)-(21d). In the frequency region where AHq = only AHq^ 
and AH^ contribute to the total linewidth, the shape of the linewidth is mainly dominated 
by AH^. This prediction is a new result. The behavior AH^ oc obtained in our model for 
high frequencies, is in contrast to conventional ferromagnetic resonance including only the 
sum of a Gilbert part linear in frequency and a two-magnon contribution which is saturated 
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at high frequencies. So far, experimentally the frequency ranges from 1 GHz to 225 GHz, 
see Let us point out that the results presented in Fig. [7] can be adjusted in such a 

manner that the Gilbert contribution will be inoperative at much higher frequencies by the 
appropriate choice of the model parameters. Due to this fact we suggest an experimental 
verification in more extended frequency ranges. Another aspect is the observation that 
excitations with a nonzero wave vector might represent one possible retardation mechanism. 
Regarding Eqs. ( 21a[ )-( 2ld) ) retardation can also influence the linewidth in case = 



(i.e. P = and Fi = Fq). Only if r = the retardation effects disappear. Therefore let us 
consider the time domain of retardation and its relation to the Gilbert damping. The Gilbert 
damping and the attenuation due to retardation can be considered as competing processes. 
So temporal feedback can cause that the Gilbert contribution disappears. In the same 
sense the Bloch contribution is a further competing damping effect. In this regard temporal 
feedback has the ability to reverse the dephasing process of spin waves based on Gilbert and 



Bloch damping. On the other hand the retardation part Ar in Eq. (21d) is always positive 
for r > 0. Thus, the retardation itself leads to linewidth broadening in ferromagnetic 
resonance and consequently to spin damping. Whether the magnitude of retardation is able 
to exceed the Gilbert damping depends strongly on the frequency. With other words, the 
frequency of the magnetic excitation 'decides' to which damping mechanisms the excitation 
energy is transferred. Our calculation suggests that for sufficient high frequencies retardation 
effects dominate the intrinsic damping behavior. Thus the orientation and the value of the 
magnetization within the retardation time r plays a major role for the total damping. 

Generally, experimental data should be fit according to the frequency dependence of the 
linewidth in terms of Eqs. (20)- (22). To underline this statement we present Fig. [s] In this 
graph we reproduce some results presented in [7J for the case Bh = ©m = 0. To be more 
specific, we have used Eq. (94) in [7j which accounts for the two-magnon scattering and 
the parameters given there. As result we find a copy of Fig. 4 in [7] except of the factor 
2 / v^- Further, we have summed up the conventional Gilbert linewidth oc / with the Gilbert 
damping parameter ai = 0.003. This superposition yields to the dotted line in Fig. [8j The 
result is compared with the total linewidth resulting from our retardation model plotted as 
solid line. To obtain the depicted shape we set the Gilbert damping parameter according 
to the retardation model a2 = 0.0075, i.e. to get a similar behavior in the same order of 
magnitude of AiJx within both approaches we have to assume that 0:2 is more than twice 
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frequency / [GHz] 

FIG. 8. (Color online) Comparison with the two-magnon model. Frequency dependence of the total 
peak-to-peak linewidth AHt for Oh = ©m = 0, /3 = 0.5, ai = 0.003, Q2 = 0.0075, T2 = 5 x lO^^s, 
r = 1.22 X W~^^ s and Tq = 1.2. Parameters taken from [7J: AttMs = 21000 G, Hs = -15000 G and 
from j37j: 7 = 0.018 GHz/G (derived from g = 2.09 for bulk Fe). The dotted line is a superposition 
of Fig. 4 in [3 reflecting the two-magnon contribution and the Gilbert contribution (denoted as 
ai in the text) linear in the frequency. 

as large compared to ai. 

Finally we discuss briefly the Gh -dependence of the linewidth which is shown in Fig. |9| 
In the upper part of the figure one observes that AH'y{Qi{) exhibits a maximum which is 
shifted towards lower field angles as well as less pronounced for increasing frequencies. The 
lower part of Fig. |9} referring to / = 10 GHz, displays that the main contribution to the total 
linewidth arises from the Gilbert part AHq. This result for / = 10 GHz is in accordance 
with the results discussed previously, compare Fig. [7j For higher frequencies the retardation 
contribution AH^ may exceed the Gilbert part. 

V. CONCLUSIONS 

A detailed study of spatiotemporal feedback effects and intrinsic damping terms offers 
that both mechanisms become relevant in ferromagnetic resonance. Due to the superposi- 
tion of both effects it results a nonlinear dependence of the total linewidth on the frequency 
which is in accordance with experiments. In getting the results the conventional model in- 
cluding Landau-Lifshitz-Gilbert damping is extended by considering additional spatial and 
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FIG. 9. (Color online) Angular dependence of the total peak-to-peak linewidtli AHj^ for various 
frequencies (top graph) and all contributions AHn for / = 10 GHz (bottom graph) with f3 = 0.5, 
a = 0.01, T2 = 5 X 10^® s, T = 1.7 x 10^^^ s and Fq = 1.1. The parameters are taken from 
[16]: IvrMs = 16980 G, Hs = -3400 G and 7 = 0.019 GHz/G. 



temporal retardation and non-conserved Bloch damping terms. Our analytical approach 
enables us to derive explicit expressions for the resonance condition and the peak-to-peak 
linewidth. We were able to link our results to such ones well-known from the literature. 
The resonance condition is affected by the feedback strength Tq. The spin wave damping is 
likewise influenced by Fq but moreover by the characteristic memory time r and the retar- 
dation length ^. As expected the retardation gives rise to an additional damping process. 
Furthermore, the complete linewidth offers a nonlinear dependence on the frequency which 
is also triggered by the Gilbert damping. From here we conclude that for sufficient high 
frequencies the linewidth is dominated by retardation effects. Generally, the contribution of 
the different damping mechanisms to the linewidth is comprised of well separated rates which 



are presented in Eqs. (20)- (22). Since each contribution to the linewidth is characterized 
by adjustable parameters it would be very useful to verify our predictions experimentally. 
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Notice that the contributions to the hnewidth in Eqs. (20)-(22) depend on the shape of 
the retardation kernel which is therefore reasonable not only for the theoretical approach 
but for the experimental verification, too. One cannot exclude that other mechanisms as 
more-magnon scattering effects, nonlinear interactions, spin-lattice coupling etc. are likewise 
relevant. Otherwise, we hope that our work stimulates further experimental investigations 
in ferromagnetic resonance. 
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